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ABSTRACT 

Recent observations have revealed that velocity dispersions of "ultra-compact 
dwarf" (UCD) galaxies are significantly smaller than those of other galaxy popu- 
lations in the Fornax and the Virgo clusters of galaxies. In order to understand the 
origin of the observed lower velocity dispersions of UCDs, we numerically investigate 
line-of-sight velocity dispersion (<ti os ) of galaxy populations with variously different 
orbits in clusters of galaxies with the total masses of M c \. We particularly investigate 
radial velocity dispersion profiles {o\ os (R)) and velocity dispersions within the central 
200 kpc of a cluster model (er m ) for galaxies with different pericenter distances (r p ) 
and orbital eccentricities (e) in the model with M c \ = 7.0 x 10 13 M Q reasonable for the 
Fornax cluster. We find that a\ os (R) and a m of galaxies with smaller r p are steeper 
and smaller, respectively, for a given initial e distribution of galaxies. For example, we 
find that a m is ~ 260 km s _1 for galaxies with r p < 50 kpc and ~ 336 km s _1 for all 
galaxies in the model with the mean e of 0.6. These results imply that the observed 
lower velocity dispersion of UCD population is consistent with the UCDs having sig- 
nificantly smaller r p than other galaxy populations in the Fornax. We discuss these 
results in the context of the "galaxy threshing" scenario in which UCDs originate 
from nuclei of nucleated dwarf galaxies. We suggest that the observed differences in 
kinematical properties between UCDs and other dwarf galaxy populations in clusters 
of galaxies can be understood in terms of the differences in orbital properties between 
UCDs and the dwarf populations. 

Key words: globular clusters: general - galaxies: star clusters - galaxies:stellar 
content - galaxies: formation - galaxies: interactions 
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1 INTRODUCTION 

Since very compact and luminous stellar systems were dis- 
covered in the central region of the Fornax cluster of galax- 
ies (e.g, Hilker et al. 1999; Drinkwater et al. 2000a, b), 
physical properties of these systems - now referred to as 
"ultra-compact dwarf" (UCD) galaxies - have been exten- 
sively investigated by observational studies (e.g., Phillipps 
et al. 2001; Mieske et al. 2002, 2004, 2006; Drinkwater et 
al. 2003; Hasegan et al. 2005; Karick et al. 2006; Firth et al. 
2006). These observations have reported very unique proper- 
ties of UCDs, such as very compact sizes (< 100 pc), possibly 
higher mass-to-light- ratio (M/L) indicative of the presence 
of dark matter, and scaling relations of dynamical proper- 
ties (e.g., internal velocity dispersions) different from those 
of GCs (e.g., Drinkwater et al. 2003; Hasegan et al. 2005). 
Physical properties of UCDs in different clusters and groups 
of galaxies are now being investigated (e.g., Jones et al. 2006; 
Kilborn et al. 2005). 



In spite of these observational progresses, it remains un- 
clear how UCDs formed and evolved in the central regions 
of clusters of galaxies (e.g., Bekki et al. 2001, 2003a; Fell- 
hauer & Kroupa 2002). Bekki et al. (2001, 2003a) proposed 
the "galaxy threshing" scenario in which UCDs originate 
from nuclei of nucleated dwarf galaxies that had been de- 
stroyed by strong tidal fields of clusters of galaxies. Fellhauer 
& Kroupa (2002) proposed that numerous smaller young 
star clusters, such as those observed in actively star-forming 
Antennae galaxy, can finally evolve into a single UCD after 
merging of them. Although the observed structural proper- 
ties of UCDs (e.g., De Propris et al. 2005), scaling relations 
of dynamical properties of UCDs (e.g., Evstigneeva et al. 
2007), and dynamical masses including possible dark mat- 
ter halos (Hilker et al. 2007) have given some constrains on 
the above theoretical models, it has not yet been determined 
which model can be regarded as more reasonable. 

Recent spectroscopic observations have reported that 
the line-of-sight velocity dispersion (cri os ) of UCDs (cti os = 
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246 ± 25 km s _1 ) is significantly smaller than those of other 
dwarf galaxy populations (396 ± 41 km s _1 ) in the Fornax 
cluster of galaxies (Gregg et al. 2007). They also reported 
that (Ti os of UCDs is even smaller than that of GCs (ui os = 
334 ± 11 km s _1 ) around NGC 1399 derived by Dirsch et 
al. (2004). Mieske et al. (2004) have also reported that the 
bright compact objects in the Fornax cluster shows a lower 
velocity dispersion than the dwarf galaxy population. The 
lower velocity dispersion has been also observed in the six 
UCDs of the Virgo cluster of galaxies (Jones et al. 2006). 
The origin of the observed lower of ai os of UCDs in these 
two clusters is one of important unresolved problems related 
to the origin of UCDs. 

The purpose of this paper is thus to try to explain 
the observed lower velocity dispersions of UCDs based on 
somewhat idealized orbital evolution models of galaxies in 
clusters of galaxies. The galaxy threshing scenario suggested 
that nucleated dwarf galaxies with small pericenter distances 
(r p ) can become UCDs: galaxy populations with only a lim- 
ited range of orbital properties can be progenitors of UCDs 
(Bekki et al. 2001; 2003a). We therefore investigate how ra- 
dial profiles of trios (aios(R)) and velocity dispersions within 
the central 200 kpc (a m ) for galaxies in a cluster depend 
on orbital properties of the galaxies (e.g., more circular or 
eccentric orbits) and thereby find galaxy populations that 
show a lower velocity dispersion in the cluster. 

The plan of the paper is as follows: in the next sec- 
tion, we describe the orbital evolution models of galaxies in 
a cluster of galaxies. In §3, we present the results on cr m and 
cios(R) for galaxies with different orbital properties. §4, we 
discuss the origin of flattened spatial distribution of UCDs 
in the Fornax cluster of galaxies based on our previous nu- 
merical simulations. We summarise our conclusions in §5. 



2 THE MODEL 

2.1 The threshing scenario for UCD formation 

Although there could be a number of explanations for the 
origin of the observed lower velocity dispersions of UCDs in 
clusters of galaxies, we here focus exclusively on the thresh- 
ing scenario and thereby discuss the observation. The pre- 
dictions of the threshing scenario are summarized as follows: 
(i) dE,Ns with smaller pericenter distances (r p ) and higher 
orbital eccentricities (e) are more likely to be transformed 
into UCDs, (ii) dE,Ns in the Fornax cluster can be trans- 
formed into UCDs, if the pericenter distances (r p ) of their 
orbits (with respect to the center of the cluster) are smaller 
than the "threshing radius" (-Rth), (iii) Rth is typically ~ 50 
kpc for dE,Ns that have B-band magnitudes (Mb) of « — 16 
mag and are embedded in dark matter halos with large cores 
(see Figure7 in Bekki et al. 2003a), and (iv) R th depends on 
Mb for a given cluster mass in such a way that R t h is smaller 
for larger (i.e., fainter) Mb- 

The proposed threshing radius has been discussed in 
the context of the observed compact spatial distribution of 
UCDs in the Fornax cluster (Bekki et al. 2003a). The or- 
bital properties of the UCDs transformed from dE,Ns with 
f p < Rth however have not been discussed at all even in 
a qualitative way. The threshing scenario implies that or- 
bital properties might well be different between UCDs ( "de- 
stroyed populations") and the presently observed dwarfs 



( "survived ones" ) and thus result in the observed kinemati- 
cal differences. In the present paper, we consider that the ob- 
served line-of-sight velocity dispersions of UCDs and dwarf 
ellipticals (dEs) can reflect the intrinsic differences in orbital 
properties between these two dwarf galaxy populations. We 
thus investigate kinematical properties of different galaxy 
populations with different r p and e in the Fornax cluster 
model. 



2.2 Orbital evolution of galaxy populations in the 
Fornax cluster 

We investigate orbital evolution of "galaxies" represented 
by test particles in fixed gravitational potentials reasonable 
for clusters of galaxies. We adopt orbital evolution models 
similar to those used by Bekki et al. (2003a) for better un- 
derstanding the origin of the spatial distributions of UCDs 
in clusters and groups. We consider that particles with the 
total number of N orbit a cluster dominated by a massive 
dark matter halo with the total mass of M c \. To give our 
model a realistic radial density profile for the dark matter 
halo of the cluster, we base it on the predictions from the 
standard cold dark matter cosmogony (Navarro, Frenk, & 
White 1996, hereafter NFW). The NFW profile is described 



p(r 
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(r/r s )(l + r/r s y 



(1) 



where r, po, and r s are the distance from the center of the 
cluster, the central density, and the scale-length of the dark 
halo, respectively. 

We mainly investigate the "Fornax" model with M c \ — 
7.0 x 10 13 M Q , r s of 83 kpc, and "c" parameter of 12.8, be- 
cause we are interested in the observed kinematics of UCDs 
and other galaxy populations in the Fornax. The adopted set 
of parameters can be consistent with the X-ray observations 
of Jones et al. (1997) and the mass estimation of the Fornax 
cluster derived from kinematics of the cluster galaxies by 
Drinkwater et al. (2001). We do not consider gravitational 
influences of cluster member galaxies including NGC 1399 in 
the Fornax model, because the influences can be very minor 
owing to the total mass of the galaxies much smaller than 
that of the background dark matter halo. 

For the projected radial number distribution of particles 
(S(i?)) within the cluster, we adopt a King profile with a 
core radius ~ 0.6 times smaller than the cluster scale radius, 
r s , in each cluster model (Adami et al. 1998). In order to 
derive the three dimensional (3D) density field (p(r), where 
r is the distance from the center of the cluster) from £(-R), 
we can use the following formula (Binney & Tremaine 1987): 



p{r) 



7T / 

J r 
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dR y/R2 - r 2 ' 
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We numerically estimate the spherical symmetric p(r) profile 
for a given £(i?) in a model. 

Each particle is given an initial orbital eccentricity (e) 
and a pericenter distance (r p ) for a given cluster potential. 
We give each particle e so that the e-distribution for N parti- 
cles has a Gaussian with the mean e value and the dispersion 
of e being e m ,o and ao, respectively. Since we adopt spher- 
ical distributions of particles, we do not introduce initial 
anisotropy in velocity dispersion for the x, y, and z direc- 
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Table 1. Model parameters 



Model no. 


e m ,o a 


-0 b 


comments 


Ml 


0.6 


0.2 


standard model 


M2 


0.5 


0.2 


more circular orbits 


M3 


0.7 


0.2 


more eccentric orbits 


Mi 


0.6 


0.1 




M5 


0.6 


0.3 





a The initial mean of orbital eccentricities (e) for particles. 
b The initial dispersion in e for particles. 

tions: the velocity ellipsoid of particles can be anisotropic 
in radial (r) direction only. In this calculation, we consider 
that the mean orbital eccentricity of galaxies in a cluster 
(em,o) is 0.6, which is consistent with recent high-resolution 
cosmological simulations (Ghigna et al. 1998). We however 
investigate different five sets of models with 0.5 < e mj o < 0.7 
and 0.1 < cr < 0.3. 

We investigate radial profiles of (projected) line-of-sight 
velocity dispersions (ai os (R)) and velocity dispersions within 
200 kpc (am) for "selected particles" and all ones. We in- 
vestigate a m within 200 kpc, mainly because all UCDs are 
located within the central 200 kpc of the Fornax (e.g., 
Drinkwater et al. 2003; Gregg et al. 2007). We consider that 
particles with different r p and e can have different kine- 
matics, and accordingly we investigate a\ OB (R) and a m for 
"selected particles" with r p (e) being within a certain pa- 
rameter range. We introduce two key parameters; (1) r t h 
which r p of selected particles are smaller than and (2) e t h 
which e of the particles are larger than. We consider that r t h 
and e t h are key, because our previous simulations of UCD 
formation based on the threshing scenario suggest that these 
two are important for UCD formation. 

We investigate the models with 25 kpc < r t h < 175 kpc 
and 0.5 < e th < 0.7 for a given e mj o and <tq. We mainly show 
the results of the "standard" model (Ml) with e mi o = 0.6 
and ao = 0.2, because the adopted parameters are the most 
reasonable in the present study. The parameter values of 
other four models (M2, M3, M4, and M5) are given in the 
Table 1. Since most UCDs are located at 50 kpc< R <200 
kpc in the Fornax cluster (Gregg et al. 2007), we estimate a m 
for particles with 50 kpc< R <200 kpc. For conveniences, 
a m for selected and all particles are referred to as er m and 
0"m, an, respectively, in the present study. 

N is set to be 10 4 for all models so that the error bar 
in each i-th radial (R) bin due to the smaller number of 
particles (i.e., ~ \f\[Wi) can be as small as ~ 20 km s . 
For models with N as large as 10 4 , we can more clearly see 
differences in kinematics between selected and all particles 
owing to the smaller error bars, though N ~ 10 4 could be 
significantly larger than a reasonable number of galaxy pop- 
ulations and intracluster stellar systems like UCDs. 

Thus we adopt somewhat idealized orbital evolution 
models for galaxies in the cluster and thereby investigate 
kinematical differences between different orbital popula- 
tions. The results presented in this study will be useful 
and helpful not only in interpreting observational results 
on UCD kinematics but also in understanding the results of 
much more sophisticated simulations with UCD formation 
models in our future numerical studies. Orbital properties, 
ages and metallicities, spatial distributions of UCDs derived 




r P < r th ( = 50kpc) 



50 100 150 200 

Radius (kpc) 

Figure 1. Dependences of line-of-sight velocity dispersions ((Ti os ) 
on projected radii (R) for all particles (solid, red) and those with 
r p smaller than a threshold radius (r t h) of 50 kpc (dotted, blue). 
Note that ai os (R) is significantly steeper and systematically lower 
in the model with selected particles with r p < r t ^. 




50 100 150 

r th (kpc) 

Figure 2. The dependence of line-of-sight velocity dispersion of 
particles within the central 200 kpc (<r m ) on r tll (=25, 50, 75, 100, 
125, 150, and 175 km s — 1 ) in the standard model (solid). cr m of 
all particles is shown by a dotted line for comparison. 

in our cosmological simulations will be given in our future 
papers (Bekki & Yahagi 2007, in preparation). Gregg et al. 
(2007) have revealed that the spatial distribution of UCDs in 
the Fornax is quite compact and appears to be significantly 
flattened. We will also discuss the origin of the observed 
flattened spatial distribution in the context of formation of 
hosts galaxies of UCDs at very high redshifts (z > 10) based 
on the results of these future cosmological simulations. 



3 RESULTS 

3.1 The standard model 

Figure 1 clearly shows that ai os (R) is significantly steeper 
and systematically smaller in selected particles with r p < 
r t h than in all ones for the central 200 kpc of the standard 
model with r t h = 50 kpc. The steeper a\ os (R) profile in 
the particles with r p < 50 kpc is due to a large number 
of more eccentric orbits (i.e., a radially anisotropic velocity 
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Figure 3. The same as Figure 1 but for the model with r t h=150 
kpc. 



ellipsoid). The difference in a\ os (R) becomes progressively 
larger with larger R, which reflects the fact that the outer 
selected particles need to have more eccentric orbits to have 
r p < 50 kpc. The velocity dispersions within the central 
200 kpc for all particles (er m ,aii) and selected ones (a"m) are 
336 km s" 1 and 260 km s _1 , respectively, for this standard 
model. 

Figure 2 shows that cr m depends on r p < rth (25 kpc< 
rth <175 kpc) such that it is smaller for smaller rth in the 
standard model: cr m /o"m,aii is 0.54 for rth = 25 kpc and 0.99 
for r t h = 175 kpc. The physical reason for this dependence 
is that the model with smaller rth has a larger number of 
particles with more eccentric orbits in the outer part of the 
simulated region (50 kpc< R <200 kpc). Figure 3 shows that 
the standard model with rth = 150 kpc shows a very minor 
difference in cri os (R) between all and selected particles, which 
confirms that rth is a key parameter for kinematics in the 
present models for a given e mj o and <to- 

Figure 4 shows that <r m depends very weakly on e t h for 
0.5 < eth < 0.7 in the sense that <r m is smaller for larger eth 
(i.e., more eccentric orbits). It is confirmed that ai os (R) is 
steeper in the selected particles with e > e t h than in all ones 
for the central region (R < 20 kpc). These results strongly 
suggest that rth is a much more important parameter than 
eth f° r velocity dispersions in the standard model. These 
results in Figures 1—4 imply that the observed lower velocity 
dispersions of UCD populations in clusters of galaxies have 
physical origins: the lower dispersions are unlikely to result 
from some uncertainties in observational methods to derive 
(7m for a smaller number of UCDs. 
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Figure 4. The dependence of <r m on e t h (=0.5, 0.55, 0.6, 0.65, 
and 0.70) in the standard model (solid). cr m of all particles is 
shown by a dotted line for comparison. 
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Figure 5. The same as Figure 2 but for different models: M2 with 
e m,o = 0.5 and oo = 0.2 (solid, magenta), M3 with e m o = 0.7 and 
<jq = 0.2 (dotted, red), M4 with e m o = 0.6 and oo = 0.1 (dashed, 
blue), and M5 with e mi o = 0.6 and on = 0.3 (long-dashed, green). 



cally smaller in the selected particles with r < rth than in all 
ones for models with different e m ,o and ero, though the de- 
tails of the profiles depends slightly on e m .o and oq. Thirdly 
dependences of cr m and a\ os (R) on eth are much weaker than 
those on rth f° r different models, which suggests that rth is 
a more important parameter for kinematics of galaxy pop- 
ulations. These results for different models are thus all con- 
sistent with those derived for the standard model. 



3.2 Parameter dependences 

The three important parameter dependences of <r m and 
(T\ OB (R) on initial orbital properties of particles in cluster 
models (i.e., e m ,o and <to) are briefly summarized as follows. 
Firstly, the dependence of a m on rth derived in the stan- 
dard model is seen in the models with different e mj o and oo. 
Figure 5 clearly shows that irrespective of e m ,o and ctq, o"m 
is smaller for smaller rth, which confirms that rth is a key 
parameter for kinematical properties of galaxy populations 
in clusters. 

Secondly, ai os (R) is significantly steeper and systemati- 



4 DISCUSSIONS 

4.1 UCDs originate from destroyed building 
blocks of clusters ? 

The present study has first shown that galaxies with smaller 
r p (< 50 kpc) can show smaller velocity dispersions in the 
central 200 kpc of the cluster models. This result implies 
that the observed lower velocity dispersions in the Fornax 
and the Virgo can be due to UCDs having smaller r p (< 50 
kpc) in comparison with other galaxy populations in these 
clusters. The result also suggests that the observed lower 
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dispersions are consistent with the galaxy threshing scenario 
in which only galaxies with very small r p can become UCDs 
after the destruction of the galaxies via strong cluster tidal 
fields (Bekki et al. 2001; 2003a). Thus the observed lower 
velocity dispersions of UCD populations can provide a clue 
to the formation processes of UCDs in clusters. 

It should be however stressed that the observed lower 
dispersions of UCDs can not rule out the alternative sce- 
nario by Fellhauer & Kroupa (2002) in which UCDs can be 
formed by merging of numerous smaller clusters in starburst- 
ing galaxy mergers ("cluster merger scenario"): if young 
UCDs formed in galaxy mergers are tidally stripped to have 
more radial orbits and smaller r p (for some physical rea- 
sons) in the central regions of clusters, the observed lower 
dispersions would be explained equally. Since galaxy merg- 
ing is highly unlikely after the formation of clusters at low 
redshifts (e.g., Ghigna et al. 1998), galaxy merging that can 
form UCDs probably needs to occur in the very early stages 
of clusters of galaxies. Thus, although the observed lower 
velocity dispersions of UCDs are consistent with the thresh- 
ing scenario, they still appear not to be evidences strong 
enough to determine which of the above two scenarios is 
more reasonable. 

4.2 Possible kinematical differences between 
UCDs and ICGCs 

Recent observational studies of GCs in clusters of galax- 
ies have suggested that there can be a population of GCs 
that are bounded by cluster gravitational potentials rather 
than those of cluster member galaxies (e.g., West et al. 1995; 
Bassino et al. 2003; Jordan et al. 2003; Tamura et al. 2006; 
Bergond et al. 2007; Williams et al. 2007). Previous numeri- 
cal simulations showed that these intracluster GCs (ICGCs) 
can be formed from tidal stripping of GCs that are located 
in the outer parts of galaxies orbiting clusters of galaxies 
(e.g., Bekki et al. 2003b). Given that the UCDs can be the 
nuclei of destroyed (thus defunct) galaxies, spatial distribu- 
tions and kinematics of ICGCs and UCDs in a cluster can 
be significantly different. 

Bekki & Yahagi (2006) suggested that ICGCs can show 
more isotropic line-of-sight velocity dispersions, which can 
be in a contrast with the possibly (radially) anisotropic ones 
of UCDs suggested by the galaxy threshing scenario (Bekki 
et al. 2003). Although observational studies on structural 
and kinematical properties of ICGCs in nearby clusters have 
just started (e.g., Tamura et al. 2006; Bergond et al. 2007), 
these dynamical properties for relatively faint ICGCs have 
not yet been derived for the entire regions of the clusters. 
It is therefore observationally unclear whether spatial distri- 
butions and kinematics within clusters are different between 
UCDs and ICGCs. Future observations on the differences in 
aios{R) between UCDs and ICGCs may well help us to un- 
derstand the differences in formation processes (e.g., thresh- 
ing vs stripping) between these two intracluster populations. 



5 CONCLUSIONS 

We have mainly investigated orbital evolution of galaxies 
with different orbital parameters (i.e., e and r p ) in the cen- 
tral region (R < 200 kpc) of the Fornax cluster model in 



order to understand the origin of the observed lower veloc- 
ity dispersion of UCDs in the Fornax. We have introduced 
two key parameters r t h and e t h and thereby investigated cr m 
and cri oa (i?) for galaxies with r p < r th (or with e < e th ) 
for a given e mj o and <jq. We have particularly focused on the 
kinematical differences between all galaxies (e.g., cr m , a ii) and 
those selected with r p < r th (e.g., a m ). We summarise our 
principle results as follows. 

(1) <Ti OB (R) is steeper in selected galaxies with r p < rth 
than in all ones for the standard model with e mj o = 0.6 
and oo = 0.2. a m is significantly smaller than cr m>a n in the 
standard model and the differences between <r m and cr m , a ii 
depend on rth in such a way that they can be larger for 
smaller rth for 25 kpc< rth <175 kpc. 

(2) aios(R) is steeper in selected galaxies with e p > eth 
than in all ones for the standard model. Although a m of the 
selected galaxies with e p > eth is slightly smaller than cr mj aii 
in the standard model, the difference between the two can 
be quite small (~ 20 km s _1 ) for 0.5 < eth < 0.7. 

(3) The above results derived for the standard model 
are seen in the models with different e m ,o and (To, which im- 
plies that lower velocity dispersions of UCDs can be due to 
their smaller r p in comparison with other galaxy populations 
in the Fornax. Thus the observed kinematical differences in 
UCDs and other dwarf galaxy populations are due to dif- 
ferences in intrinsic orbital properties between these galaxy 
populations in the Fornax cluster. 

Although the observed significantly flattened spatial 
distribution of UCDs in the Fornax cluster can potentially 
have important implications on dwarf galaxy formation in 
the cluster (Gregg et al. 2007), we did not discuss the origin 
of the distribution. We plan to discuss the observed distribu- 
tion in the context of the picture that UCDs originate from 
galaxies embedded in dark matter halos already virialized 
at z > 10 in the proto-cluster environment (Bekki & Yahagi 
2007, in preparation). 
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